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Linear (spectro-)polarimetry is usually performed using separate photon flux 
measurements after spatial and/or temporal polarization modulation. Such 
classical polarimeters are limited in sensitivity and accuracy by systematic 
effects and noise. We describe a spectral modulation principle which is based 
on encoding the full linear polarization properties of light in its spectrum. 
Such spectral modulation is obtained with an optical train of an achromatic 
quarter-wave retarder, an athermal multiple-order retarder and a polarizer. 
The emergent spectral modulation is sinusoidal with its amplitude scaling 
with the degree of linear polarization and its phase scaling with the angle 
of linear polarization. The large advantage of this passive set-up is that 
all polarization information is, in principle, contained in a single spectral 
measurement, thereby eliminating all differential effects that potentially 
create spurious polarization signals. Since the polarization properties are 
obtained through curve-fitting, the susceptibility to noise is relatively low. 
We provide general design options for a spectral modulator and describe the 
design of a prototype modulator. Currently, the set-up in combination with 
a dedicated retrieval algorithm is able to measure linear polarization signals 
with a relative accuracy of 5%. (c) 2009 Optical Society of America 
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1 



1. Introduction 

Polarimetry is a powerful remote-sensing method, but its implementation can be technically 
challenging. Classically, the complete information on the linear polarization properties of 
light is obtained from at least three independent intensity (i.e. photon flux) measurements 
that are "filtered" for a particular polarization state. As an example, we discuss full linear 
polarimetry with a four-fold modulation scheme: 

si = t{i + Qy, 

52 = t{I-Q)- 

53 = t{I + U)- 

54 = t{I-U), (1) 

with t as the transmission factor, which is ideally identical for all measurements. The linear 
Stokes parameters Q and U and the intensity / (all multiplied by the transmission t) are 
trivially obtained from the four intensity measurements Si ,4. Usually, the normalized pa- 
rameters Q/I and U/I are determined, since those are independent of t. The polarization 
modulation ensures that for each measurement the correct polarization state is transmitted 
by a polarizer. Two different classical modulation principles are discerned: spatial and tem- 
poral modulation. In the case of spatial modulation, the light is split into four beams by, e.g., 
polarizing beamsplitters that select light with linear polarization at 0° and 90° (for determin- 
ing Stokes Q from 51,2) and ±45° (Stokes U). With temporal modulation, the measurements 
S1..4 are obtained sequentially with a single beam by, e.g., rotating the polarizer or rotating a 
half-wave plate in front of a fixed polarizer. Since the four intensity measurements are inde- 
pendent, polarimetry using spatial or temporal modulation is sensitive to various systematic 
differential effects as well as to noise. Spatial modulation requires four different optical paths 
and four detectors or four different parts of an imaging array to measure the four signals. 
Therefore, differential alignment and transmission of the four beams is very critical for a pre- 
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cise polarization measurement. Also, limited calibration of the gain or flat-field of the imaging 
array introduces errors in the polarimetry. On a more practical side, the division into three 
or four beams significantly increases an instrument's size and mechanical complexity. As for 
temporal modulation, differential effects that reduce the polarimeter's performance include 
time-variations of the source or of the instrument's pointing. For instance, a remote-sensing 
instrument that is scanning a source in time as well as modulating polarization in time is 
prone to spurious polarization signals. Also, varying atmospheric properties (such as trans- 
mission changes or "seeing" due to turbulent refractive index variations) or vibrations of the 
instrument create such spurious signals. Temporal modulation requires an active component 
such as a rotating wave plate or a liquid crystal device, which consume power and may fail 
at some point during the instrument's life-time. Many differential errors are eliminated to 
first order by combining spatial and temporal modulation with the so-called beam-exchange 
method [1]. Nevertheless, the mechanical complexity of a beam-exchange polarimeter is even 
larger than for a spatial or temporal modulation polarimeter. Moreover, all the described 
methods for polarimetry are sensitive to the random noise in the individual measurements 
■51..4, particularly when the degree of polarization to be measured is very low. The difference 
measurements Si — S2 and S3 — S4 are then usually dominated by random noise. 

Here we present a new method for measuring linear polarization as a function of wavelength 
that is much less sensitive to differential effects. If the polarization information cannot be 
stored in the spatial or temporal dimension, and if the instrument is designed to measure 
the spectral dependency of the light as well (spectropolarimetry), then it is possible to store 
the polarization information in the spectral dimension. Such a method, alternately known 
as channeled spectropolarimetry or polarimetric spectral intensity modulation, was reported 
independently by Oka & Kato (1999) [2] and lannarilli et al. (1999) [3,4], but their set-ups 
measure the complete Stokes vector including circular polarization, whereas our method only 
applies to linear polarization. It is therefore convenient to describe the polarization properties 
in terms of the degree of linear polarization (Pl) and the angle of linear polarization {(pi), 



3 



which are obtained from the following coordinate transformation: 



(2) 



I U 

<pL = -arctan— . (3) 

A "spectral" modulation for linear polarization is obtained through encoding the full linear 
polarization properties into a sinusoidal variation of the intensity spectrum such that the 
amplitude of the sinusoid scales with and its phase scales with 0^. A spectral modulator 
with such an output can be constructed from standard, passive optical components: a com- 
bination of an achromatic quarter-wave retarder, an athermal multiple-order retarder and a 
polarizer, as is elaborated in sections 2 and 3. The major advantage of this spectral modu- 
lation principle is that the measurable polarization properties are, in principle, contained in 
a single intensity spectrum recording, which eliminates sensitivity to any differential effect. 
Moreover, since Pl and (pL can be obtained through a curve-fit of the sinusoidal spectral 
modulation, the susceptibility to noise is lower compared to spatial and temporal modu- 
lation. Of course the implementation of spectral polarization modulation comes at a cost: 
the spectral resolution of the instrument needs to be much larger than the required spectral 
resolution of the data-product to carry the additional polarization information in the spec- 
trum. However, compared to the Oka & Kato method, which encompasses three sinusoidal 
modulations with different spectral periodicities, the method described here minimizes the 
increase in spectral resolution for linear spectropolarimetry. The Oka & Kato method is lim- 
ited by the thermal variation of the employed multiple-order crystals [5]. In section 3.B we 
describe a way to create "athermal" retarders. But even with residual thermal variations of 
the multiple-order retarder, the measurement of Pl (which is often the key observable) is, 
in principle, unaffected. 

In this paper, we aim to both describe the general principle of spectral polarization mod- 
ulation (section 2), as well as the technical implementation thereof. Optical design consid- 
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erations for the construction of a spectral modulator are discussed in section 3. We created 
a dedicated algorithm that successfully disentangles the polarization information from the 
intensity spectrum. The algorithm and first results from a prototype spectral modulator are 
presented in section 4. In section 5 we discuss the systematic instrumental errors pertaining 
to the retrieved values of Pl and (pL and methods to calibrate those. 

We foresee a number of applications for full linear spectropolarimetry by spectral modula- 
tion. Any remote-sensing measurement that benefits from the additional value of obtaining 
the linear polarization properties as well as the spectrum qualifies. In particular, interesting 
applications include the characterization of aerosols by analyzing scattered light [6] and the 
detection of skin cancer [7], land mines [8], camouflaged objects [9] and of disease in crops 
and vegetation [10]. In all these cases polarization is created or modified by large-angle scat- 
tering or reflection of natural or artificial light. With a spectral modulator, any spectrometer 
(spectrograph or hyperspectral imager) can be transformed into a spectropolarimeter, albeit 
at reduced spectral resolution. 

2. Spectral modulator principle 

The spectral modulator for full linear polarimetry creates a sinusoidal variation of the spec- 
trum for which the amplitude scales with the degree of linear polarization {Pl) and the phase 
is determined by the angle of linear polarization Such a spectral modulation can be 

obtained by a combination of an achromatic quarter-wave retarder, an athermal multiple- 
order retarder and a polarizer. The schematic set-up including the mutual orientations of the 
individual components is shown in Fig. l.a. The quarter- wave retarder's fast axis and the 
polarizer are oriented in parallel, although a 90° rotation of both components does not af- 
fect the modulation principle (see below). The athermal multiple-order retarder has its axes 
oriented at ±45° from the other two components. It consists of two crossed wave plates from 
different crystal materials to make it athermal as explained in section 3.B. The retardance 
of the combined multiple-order crystals is described by 5(A, T) (see Eq. 11). With the set-up 
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of Fig. l.a, any spectrum So(A) is spectrally modulated as: 
SW = IsoiX) ■ 

with the ± sign depending on the orientation of the polarizer. The orientation of (pL at the 
source depends on the ±90° orientation of both retarders. Note that the modulation is not 
strictly periodic in wavelength, but in 1/A (oc frequency). The optimum value of 6 depends 
on the required spectral resolution. A typical value for 6 for a spectral modulation period 
of ~ 20 nm (which roughly corresponds to the final data-product spectral resolution) in the 
visible is 25 /im, as obtained after optimization using the retrieval algorithm introduced in 
section 4. The spectrum So(A) is defined as: 

soiX) = /o(A)t(A), (5) 

where Iq is the intensity spectrum of the light that enters the spectral modulator. t{X) 
includes the transmission of the instrument, the detector gain and the analog-to-digital 
conversion. 

The validity of Eq. 4 is easily verified with the Mueller matrix calculus. As an illustration. 
Fig. l.b shows the consecutive operations of the polarization components on the Poincare 
sphere. The measurement phase space of linear polarization is spanned by the {Q/I,U/I) 
plane on which fully polarized light traces out the equator of the sphere and completely un- 
polarized light is represented by the center of the sphere. The (achromatic) quarter-wave re- 
tarder has its axes in the direction of ±Q and rotates the {Q/I,U/I) plane to the {Q / 1, V/ 1) 
plane (with Stokes V representing circular polarization). The multiple-order retarder has its 
axes in the ±f/ direction and rotates all points in the {Q/I, V/I) plane around the U/I axis. 
Since multiple-order crystal retarders are by nature very chromatic, the amount of rotation 
depends much on wavelength. Finally, the polarizer in the ±Q direction makes a projec- 
tion of all points in the sphere onto the Q/I axis, the position on which corresponds to the 



1±Pl(A)-cos 



27r-^(A,T) 



2 ■ 0l(A) 



(4) 
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relative amount of light (from to 1) that is transmitted by the polarizer. The action of 
the multiple-order retarder is thus transformed into the sinusoidal intensity spectrum. This 
is essentially the same principle as employed in Lyot filters [11]. From this illustration it 
is apparent that the degree of linear polarization Pl directly scales the amplitude of the 
sinusoidal modulation: fully polarized light creates spectra with maximum modulation am- 
plitude < ///q < 1, whereas completely unpolarized yields a fiat spectrum of ///q = \. 
A variation of the angle of linear polarization 0l is represented by a rotation around the 
V/ 1 axis of the Poincare sphere. This rotation is maintained by the action of the achromatic 
quarter-wave retarder and translated into a phase shift of the sinusoidal modulation by the 
other two components. Note that a similar spectral modulation can, in principle, be achieved 
by replacing the two (linear) retarders by an optically active material (circular retarder). But 
in order to create a narrow spectral periodicity of ~ 20 nm it would require unreasonably 
thick optically active crystals. 

The polarizer of the spectral modulator can also be a polarizing beam-splitter. Eq. 4 implies 
that the sum of the spectra emerging from a polarizing beam-splitter is ideally equal to the 
intensity spectrum at full resolution (!), provided that the transmission ratio ti(A)/t2(A) (see 
Eq. 5) between the two beams is known. Dividing both measured spectra by the intensity 
spectrum then yields two independent modulation signals scaled as < s/sq < 1. This 
represents a straightforward method to disentangle the polarization information from the 
intensity spectrum. Particularly in cases where narrow features of the intensity spectrum mix 
up with the polarization modulation, such a dual beam set-up is recommended. However, if 
the intensity spectrum is sufficiently smooth, an algorithm can be devised that accurately 
retrieves the intensity spectrum as well as Pl and 0l as described in section 4, and a single 
beam set-up suffices. 

As with any polarimeter, the modulator is best located as early in the beam as possible 
in order to minimize the number of optical components that modify the polarization of the 
source under investigation and/or introduce instrumental polarization. Ideally, the spectral 
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modulator is located in the entrance pupil of a spectropolarimetric instrument. The objective 
lens(es) are then positioned after the polarizer or polarizing beam-splitter to image the source 
onto the entrance (slit) of the spectrometer. 

3. Technical implementation 

The technical implementation for the spectral polarization modulator as schematically de- 
picted in Fig. l.a depends much on the scientific requirements of a particular instrument. 
Different options exist for the optical components based on the required wavelength range, 
the field-of-view (FOV) and the polarimetric accuracy. Also, mechanical issues can impose 
constraints on the optical materials and components (e.g. robustness, clear aperture, com- 
plexity). In this section we discuss different design options for the optical components of a 
spectral polarization modulator. We also present our chosen design for a spectral modula- 
tor that operates in the extended visible range in combination with a 1 nm resolution slit 
spectrograph. 

3. A. Achromatic quarter-wave retarder 

Quarter-wave plates based on birefringent crystals are inherently chromatic. So-called 
(super-) achromatic wave plates can be constructed by using pairs of different crystals [12] or 
a stack of wave plates at different angles [13] or a combination of both these methods. Nev- 
ertheless, residual wavelength variations of the retardance and of the fast axis angle persist 
for these wave plates. 

A second variety of retarders is based on total internal reflections (TIRs) inside glass 
rhombs for which the retardance due to a single TIR is given by [14]: 



The major advantage of TIR retarders is the achromaticity of the retardance, which only 
varies with the wavelength variation of the refractive index n{X). The most famous example 
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of a quarter-wave retarder based on TIRs is the Fresnel rhomb, which exhibits ±2° retardance 
variation over the visible range when made from BK7 glass. Unfortunately, Fresnel rhombs 
for the visible range cannot be manufactured out of fused silica, which is a favorable material 
because of its low amount of birefringence induced by internal stresses or by temperature 
variations [14]. It is possible to tweak the overall retardance of a Fresnel rhomb by means of 
a coating on the TIR surface [15], but temperature variations can still lead to unacceptable 
variations of the rhomb's retardance. A K-prism based on 3 TIRs can be made out of fused 
silica and has the additional advantage of the lack of lateral beam shift, but it yields a much 
thicker optical component than a Fresnel rhomb. The FOV behavior of both a Fresnel rhomb 
and a K-prism is very anamorphic: the retardance changes by as much as 9 x 10~^ radians 
for a 1° change in beam angle in the plane of the TIRs whereas it is almost constant for 
beam angle changes in the perpendicular direction. This behavior does however not severely 
limit the performance of a spectral modulator in combination with a long-slit spectrograph 
in the case that the slit is aligned with the "good" direction of the Fresnel rhomb or K-prism. 
Another type of quarter-wave TIR prism from fused silica can be designed if the output beam 
direction is allowed to deviate from the input direction as sketched in Fig. 2. This way, also 
the FOV behavior of the prism can be optimized. For a retroreflecting prism, the equations 
for determining the three TIR angles (a,/?, 7) are: 



< a,/5,7 < 90°, 



(7) 



with 6c the critical angle for TIR; 



a + /3 + 7 = 180°; 
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(5TiR(a, Ao) + 5tir(/5, Ao) + 5tir(7, Aq) 



TT 



(9) 



2 



d^TiR(«, Aq) _ d(5TiR(/3, Aq) d(5TiR(7, Ap) 
d^ d^ d^ 



0. 



(10) 
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The minus sign in the last equation is due to the mirroring effect from the first TIR {a; see 
Fig. 2). A solution to Equations 7-10 is found for the case of fused silica with a = 82.5°, 
(3 = 50.0°, 7 = 47.5° at Aq = 500 nm. Its optimized FOV yields a retardance variation of 
—3 X 10~^ radian for a ±2° variation of the beam angle in the plane of the rhomb. Note that 
variation of the incidence angle causes significant beam shift as indicated in Fig. 2. 

3.B. Athermal multiple- order retarder 

The retardance of a crystal wave plate is usually very dependent on temperature, particularly 
in case of multiple-order retardance. Inspection of Eq. 4 shows that such a temperature 
variation of (5(A, T) has an identical effect on the spectral modulation as a variation of (pL (i-e. 
rotation of the linear polarization direction) has. This thermal dependency of the multiple- 
order retarder therefore must be eliminated to unambiguously measure 0l. Active thermal 
control is an option, although it increases the modulator's (thermo-)mechanical complexity 
and power consumption. A second option is to create an "athermal" wave plate by combining 
two (or more) plates from different crystals with different thermo-optic constants, such that 
their thermal variations cancel (and their retardances don't). Note that the measurement of 
Pl is to first order unaffected by (residual) thermal variations of 5(A,T). 

3.B.I. Design options 

To create an "athermal" multiple-order retarder, two different materials need to be found 
such that for a certain thickness ratio the residual thermal dependency of the combined 
retardance is minimal over the required wavelength range. The effective retardance of a 
single retarder plate is a strong function of incidence angle [14]. Preferably, for the crystal 
plate combination the FOV behavior is optimized. Combinations of three or more crystal 
plates can be considered in order to gain additional degrees of freedom in the design space. In 
the case of two plates i = 1, 2 with thicknesses di and birefringence ne,i — rio^i, the combined 
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Fig. 1. (a) Schematic set-up of the spectral modulator. The solid and dashed 
lines of the retarders represent the (orthogonal) fast and slow axis respectively, 
(b) Illustration of the spectral modulation principle on the Poincare sphere. 
Snik_AO105143_figl.eps 




Fig. 2. Design of a retroreflecting quarter- wave retarder prism based on 3 TIRs 
in a fused silica rhomb. The path of light with an angle of 2° from normal 
incidence is indicated in light grey. Snik_AO105143jig2.eps 
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retardance S{X,T) is given by: 



6{X,T)=6i{X,T)±62{\,T)= ne,i(A, T) -r2o,i(A, T) ■ ± n^.slA, T) -no,2(A, T) -d^, (11) 



with the ± sign for an adding/subtracting combination of crystals. The absolute value of the 
birefringence is taken in order to describe the wave plates' alignment between their fast axes 
instead of their optic axes. The requirement for combined athermal behavior is formulated 
as [16,17]: 

7i-5i(Ao)±72-(52(Ao) =0, (12) 



with 



7i 



1 d5,(Ao) 1 dd. 



d(ne,j(Ao) - rio^iiXo)) 



(13) 



6i{Xo) dT didT ne,i(Ao) - r2o,i(Ao) dT 
The vast majority of common crystals exhibit negative values of 7^ [16, 18] and therefore a 
subtractive combination of crystals has to be chosen (minus sign of ± in Equations 11 and 
12) to build an athermal retarder. 

The FOV behavior of the crystal combination depends on the signs of rie^i — Uo^i and 
^£,2 — '^0,2- Assuming that the first crystal has positive birefringence (e.g. quartz or MgF2), 
the combined retardance for a subtractive combination as a function of a beam's inclination 
C, and azimuth 77 is given by [19]: 



5(A,C,^) 



5,{X) 
-S2{X) 
-52(A) 



/ cos^ Tj siv? r]^ 



n 



0,1 



1 + 





^sin^ r] 


cos^ 


1no.2 


\ ^e,2 


no,2 j _ 




^COS^ T] 


sin^ Tj \ 




\ '^e,2 


no,2 ) _ 



ne,2 - no,2 > 0; 

ne,2 - no,2 < 0. (14) 



The dependence of the retardance on the azimuth 77 can be made to (almost) disappear by 
an appropriate subtractive combination of crystals having both positive (or both negative) 
birefringence. This effectively transforms the hyperbolical curves of equal retardance on the 
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FOV into elliptical curves, which, on average, have a flatter behavior with the inclination. 
However, the combination of crystals with opposite birefringence yields a larger FOV with 
relatively constant retardance for an anamorphic slit aligned ±45° from the crystals' axes 
(i.e. in the x ot y direction in Fig. l.a) up to field angles larger than ±10°, whereas the 
acceptance cone half-angle is about 2°. Since the multiple-order retarder has to be aligned 
±45° from the axes of the quarter- wave retarder (see Fig. l.a), a combination of a Fresnel 
rhomb or K-prism and the athermal multiple-order retarder consisting of opposite crystals 
with a slit spectrograph yields the most constant behavior along the slit. Note that, in 
principle, a constant value of 6 needs to be only guaranteed over the FOV that one spatial 
pixel within the spectrometer spans. The response of the spectral modulator then needs to 
be calibrated for each spatial pixel independently. 

The full Stokes spectral modulation method of Oka & Kato [2] would also significantly 
benefit from the use of athermal retarders. 

3.B.2. Thermal tests 

Using tabulated values for 7j from [18], we searched for an optimal combination of crystal 
that yields athermal behavior in the (extended) visible range from 350 to 800 nm. We selected 
the common, robust crystals quartz, MgF2 and AI2O3 (sapphire) which all have high trans- 
missions over a very large wavelength range (UV-NIR) and investigated the combinations of 
opposite crystals: quartz-sapphire and MgF2-sapphire. The resulting thermal residuals over 
the wavelength range as a function of the thickness ratio di : ^2 were found to be minimal 
for a MgF2-sapphire combination with di: 62 = 2.7. 

However, the thermo-optic constants from the literature are too inaccurate for the exact 
determination of the thickness ratio for the athermal combination. We therefore performed 
thermal tests on a combination of MgF2 and sapphire crystals with a variable thickness ratio. 
A 1.1 mm thick sapphire crystal was mounted in a Soleil-Babinet mount from B. Halle in 
(subtractive) combination with a wedged MgF2 crystal with a thickness varying from 1.2 to 
4.6 mm. This set-up was positioned with the crystal axes at ±45° between crossed Glan-laser 
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polarizers and put in a thermal chamber. Light from a halogen source outside the chamber 
was fed into the set-up in the chamber by means of an optical fiber. A microscope objective 
and two diaphragms were used to produce a collimated beam of ~1 mm diameter through 
the polarizers and the crystals, which is small enough to assume a constant thickness ratio 
over the beam. Finally, another microscope objective feeds the light into a second fiber that 
is connected to an Ocean Optics USB2000 spectrograph (also outside the chamber) with a 
spectral range of 180-880 nm with 0.34 nm/pixel and a spectral resolution of 0.6 nm. The 
temperature inside the chamber was varied from 5 to 50°C in three equal steps for a number 
of different thickness ratios. 

This set-up is essentially the same as depicted in Fig. l.a. The first polarizer that de- 
termines the input polarization to the modulator and its orientation is fixed and oriented 
in the same direction as the quarter-wave retarder would be. This input polarization is an 
eigenvector of the quarter-wave retarder, which is therefore not implemented in this set-up. 
Because 0l is fixed, a measured variation of that angle is directly related to the temperature 
variation of the multiple-order retarder (see Eq. 4): 

The angle of linear polarization (pi is retrieved by the algorithm presented in section 4. For 
a range of thickness ratios between MgF2 and sapphire and for wavelengths from 550 to 750 
nm (where signal/noise was optimal), the variation of 6 with temperature was obtained after 
linear fits of the data to Eq. 15. The results, after averaging over 30 nm wavelength bins, are 
plotted in Fig. 3. a. It is clear that the data cross the zero of the vertical axis which represents 
the thermal variation of the combined retardance. The exact values for the thickness ratio 
for which the data cross zero are determined from second order polynomial fits to the data 
and are represented by the crosses. 

For comparison, the theoretical curve based on thermo-optic constants from [18] is plotted 
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(a) 




Fig. 3. (a) Thermal test results for the combination of MgF2 (with variable 
thickness) and sapphire as a function of wavelength (550-750 nm). Various 
thickness ratios (MgF2:Al203) from 1.7 to 3.7 have been measured. The vertical 
axis represents the relative variation of the overall retardance of the crystal 
combination with temperature. The crosses represent the determined zero- 
crossings of the data for a certain wavelength. The dashed line represents the 
thickness ratio 1:2.2. (b) The corresponding theoretical curve with literature 
values for the thermo-optic constants. The dashed line here represents the 
zero-crossing of the curve. Snik_AO105143_fig3a.eps; Snik_AO105143_fig3b.eps 
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in Fig. 3.b. The general shape of the data is reproduced, but the thickness ratio values 
for which athermal behavior is obtained differ by about 0.3. The exact error bars on these 
results are difficult to estimate, but they are not very relevant to us as explained below. 
Assuming that the inaccuracy of the literature values are the main cause for the discrepancy 
between the data and the theoretical curve, we adopt the results from the measurement 
data as the new best thickness ratio di : d2 = 2.4. It is clear from both the data as well 
as the theoretical curve that exact athermal behavior cannot be achieved for the entire 
wavelength range simultaneously. The obtained thickness ratio therefore merely represents 
an optimum value for the wavelength range (350-800 nm). At the edges of the wavelength 
range the residuals are less than 1.5x10"^ radians/°C. The situation can change slightly if 
the discrepancy between Fig. 3. a and b is due to (hard-to-judge) systematic errors in the 
determined thickness ratio. In any case, the measurements prove that the set-up is perfectly 
athermal somewhere within the wavelength range. The impact of thermal residuals on the 
modulator performance and a method to mitigate these effects are discussed in section 5. 

3. C. Polarizer 

Many options exist for the choice of the polarizer or polarizing beam-splitter, mostly de- 
pending on the wavelength range. The extinction ratio should be as large as possible, hence 
crystal polarizers such as Glan or Wollaston prisms are preferred. If the spectral modulator 
is located in the entrance pupil of the instrument and is based on a dual beam system, a 
Wollaston followed by one or two objective lenses forms two separated images of the sources. 
If the spectrometer is a slit spectrograph, the separation by the Wollaston prism needs to 
be parallel to the slit direction. 

4. Retrieval algorithm 

Since the spectral modulator entangles the linear polarization information in the intensity 
spectrum, a dedicated algorithm was developed to retrieve all relevant information. We 
assume that and (pL can vary randomly with wavelength. Together with the facts that the 
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modulation is only periodic with the inverse wavelength and that S varies with wavelength, 
this seriously complicates retrieval methods based on Fourier or wavelet analyses or on a 
Hilbert transform. Instead, we fit Eq. 4 to a range of narrow wavelength windows. The size 
of these windows should be larger than one modulation period and roughly determines the 
final spectral resolution. The algorithm is slightly different for a single beam and a dual beam 
system (i.e. without or with a polarizing beam-splitter), because the intensity spectrum is 
trivially obtained in the second case by adding the two beams. We choose not to assume any 
knowledge about the physical parameters of the spectral modulator and the spectrometer, 
because the error propagation from those parameters is sometimes difficult to judge. All non- 
ideal effects are to be taken into account by calibration (see section 5). A first estimate for 
the dispersion of the retardance 6{X) is supplied to the retrieval algorithm, but this estimate 
can be improved after calibration. 
The basic steps of the algorithm are: 

1. Determination of the spectral windows. For each wavelength value measured by the 
spectrometer, a spectral window equal to the local modulation periodicity AA is de- 
termined as follows for a certain order k of the multiple-order retarder: 



The modulation periodicity therefore has a quadratic behavior for 6 ^ X. At the 
edges of the wavelength range, no windows can be created that completely contain 
measurement values, so some data is lost there, especially on the red side of the spectral 
range where the modulation period is much longer. 

2. For the single beam system the intensity spectrum at the center of each spectral window 
is estimated as the average value of the data within the spectral window. For the dual 



S = k\^{k- 1)(A + ^) ^ (A; + i)(A - ^); 



(16) 



AA 



5(1 + AV4(52)' 



(17) 
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beam system the intensity spectrum is obtained by adding both modulated spectra. 

3. The measured spectrum is divided by the obtained intensity spectrum. 

4. A first estimate for Pl is obtained by determining the minimum and maximum value 
of the normalized data within each window. 

5. An artificial reference signal is created with 0l = and Pl = 1. 

6. A first estimate for (^l in each window is obtained after cross- correlation of the nor- 
malized signal and the reference signal and a cos^(A + A<^) fit of the cross-correlation 
signal. The determined phase shift A,^ between the data and the reference signal is 
transformed into an estimate of the angle of linear polarization: 



^L,est. = A A / (18) 



AA(A) 



7. With the initial estimates for Pl and 0l, full curve-fitting of the normalized data 
within each window is performed (cf. Eq. 4). Because the solution can run away with 

we apply periodic constraints such that < </)l < vr and < < 1. The values 
of 5(A) are calculated from the known crystal thicknesses and literature values for the 
birefringence [18]. 

8. The final results for -Pl(A) and </'l(A) are smoothed with the window widths of Eq. 17 
because information at higher spectral resolution cannot be determined reliably. 

Experimental results obtained with a prototype spectral modulator consisting of a fused 
silica K-prism, a subtractive combination of 2.53 mm of MgF2 and 1.10 mm of sapphire, 
and a rotatable Glan-laser polarizer in combination with the fiber-fed 0.6 nm resolution 
spectrograph mentioned above are presented in Fig. 4 as well as the retrieved values of 
Pl(A) and </'l(A). The polarizer was used to mimic a dual beam system by rotating it back 
and forth by 90°. The retrieved values using the single beam method (dotted curves in the 



lower panel) are very similar to those of the dual beam method (solid curves). The input 
polarization was created with another rotatable Glan-laser polarizer, and therefore Pl should 
be equal to 1 and 0^ should be constant. A rotation of the input polarizer leads to the correct 
increase of the retrieved 0l(A). However, it is obvious that some deviations are observed in 
the measured values of both -Pl(A) and (pii^)-, and therefore these measurements represent 
a calibration of non-ideal effects (see next section). 

Using simulated spectral data in combination with the retrieval algorithm, we determined 
the optimum value of 5 for a required final spectral resolution of 20 nm using a spectrom- 
eter with 2 nm resolution from 350-800 nm. The minimum value for the squared difference 
between simulated and retrieved data was found for 5 ~ 25/im. The accuracy with which 
the simulated data was retrieved was within 5% relative for Pl and within ±10° for (pi 
(see [6]). It is difficult to judge what the ultimate sensitivity (i.e. the noise level of the P^ 
measurement) and accuracy of the spectral modulator can be, since this is likely dominated 
by unforeseen instrumental effects. Further measurements with (and calibrations of) the 
spectral modulator prototype are required to investigate this. 

5. Error sources and calibration 

A number of instrumental error sources can be identified a priori. We present a list of 
potential major error sources and ways to mitigate those. Fortunately the measurement 
principle is quite independent for P^ and (pL so that most errors influence only one of these 
measurements. 

5. A. Static errors 

Many alignment errors and static, non-ideal properties of the components of the spectral 
modulator can lead to a reduced efficiency of the Pl measurements. These efficiencies can 
depend on 0^ as shown in the lower panel of Fig. 4. They are also likely to vary across the 
FOV. All these effects are assumed to be stable in time and are calibrated for each pixel in 
the FOV by a measurement like in Fig. 4 as elaborated in section 5.B. A number of effects 
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Fig. 4. The upper panel shows a modulated spectrum obtained with the pro- 
totype described in the text (solid curve). The dotted curve represents the 
unmodulated intensity spectrum obtained after adding the signals from two 
orthogonal orientations of the polarizer. The normalized modulation signal is 
shown below. The lower panel shows the results for -Pl(A) and (/'l(A) as ob- 
tained by the retrieval algorithm. Offsets from the ideal values of -Pl(A) (=1) 
and 01, (A) (constant) are clearly observed. The solid curve corresponds to the 
signal shown in the upper panel using the dual beam method, whereas the 
dotted curve shows the results from the single beam method {(pL is the same 
in both cases). The dashed curves represent rotations of the input polarization 
in steps of 45°, which are reproduced by the (pi measurements. The dashed 
curve for -Pl(A) represents the measurement with the input polarization direc- 
tion coinciding with the K-prism axis, for which the measurement efficiency is 
significantly higher than for the situation with input at 45° from that (solid 
curves). Snik_AO105143_fig4.eps 
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can influence the Pl measurement: 

• Opto-mechanical alignment errors can influence the measurement efficiency in many 
ways. Individual alignment sensitivities can be assessed through the Mueller matrix 
calculus. 

• An imperfect quarter-wave retarder results in an incomplete translation of the 
{Q/I,U/I) plane to the {Q/I,V/I) plane (see Fig. l.b). The modulation amplitude 
decreases, and therefore linear polarization in the ±U direction is measured with less 
efficiency. This effect is clearly visible in the data shown in Fig. 4. 

• Pleochroism (partial polarizibility) in (one of) the multiple-order retarder crystals re- 
duces the modulation amplitude, but this depends on 0^. An extinction ratio E of the 
combined multiple-order retarder (which should be close to 1) leads to a maximum 
reduction of the measured Pl of 2\/E/{1 + E). 

• An imperfect polarizer also directly reduces the modulation amplitude. 

• The modulation amplitude is also reduced by convolution of the spectra to be measured 
with the spectrograph profile. This effect is more severe in the blue side of the spectrum 
as is seen in Fig. 4. 

• Variation of 5 over a spatial pixel or over the beam footprint on the multiple-order 
retarder leads to a smearing of modulations with slightly different periodicities, which 
leads to an effective reduction of the modulation amphtude. 

The zero point of the (pii^) measurement can exhibit offsets due to the following causes: 

• Limited opto-mechanical alignment accuracy. 

• The actual value of 6{X) is not perfectly well known. This influences the retrieval of 
0(A) as shown in the lower panel of Fig. 4. Assuming that 6{X) is constant in time and 
that the same retrieval algorithm is used for both the calibration and for the actual 
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data, the calibration determines the zero point of 0^;^ as a function of wavelength. It is 
expected that this effect also influences the retrieval efficiency of Pj,, but this is then 
consistently calibrated in the same manner. The retrieved trend of could be used 

to establish the exact values of 5(A), such that becomes fiat. 

5.B. Calibration 

The static errors listed above can all be calibrated by controlling the input polarization 
(preferably the degree and the angle of linear polarization) and analyzing the output of the 
spectral modulator. The most straightforward calibration routine consists of positioning a 
rotating (perfect) polarizer in front of the spectral modulator. The retrieved values of -Pl(A) 
and 0l(A) represent the reference data for any other observation with the spectropolarimeter. 
Since the degree of linear polarization measurement potentially depends on the angle, the 
efficiency of -Pl(A) needs to be mapped as a function of 0l(A). Care has to be taken that the 
calibration beam does not have a Stokes V component, since an imperfect spectral modulator 
may be sensitive to circular polarization. 

The calibration set-up with a rotating polarizer is not able to vary the degree of linear 
polarization, which is desirable since many spectropolarimeters are designed to measure 
signals with <C 1. Several (non-linear) instrumental effects that hamper the accurate 
measurement of low degrees of polarization may not be properly calibrated with an input 
Pl = 1. Also, the sensitivity of the polarimeter can thus be established and perhaps improved 
by identifying limiting effects. A calibration set-up which allows for full control of 1 < Pl < 
and (pL can be constructed using a 'theta cell' [20-22]. In any case, the accuracy of the 
calibration of a spectral modulator should be at least 5% of Pl and ±10° in 0^, because the 
noise levels in the retrieved values of Pl and 0l are easily suppressed below those values. 

5.C. Dynamic errors 

A number of dynamic effects can modify the performance of a spectral modulator in time. 
The influence of these effects could be detectable in the data: 
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• Bias drift and variable stray light with a measured total amount lead to a reduction 
of Pl by so/(so + 2sd). It is therefore crucial to use a detector with negligible dark drift 
and to deal with stray light by adequate baffling. 

• The transmission ratio ti(A)/t2(A) between the two beams can slightly vary over time 
e.g. due to aging of optical coatings and/or the collection of dust. This effect could 
be detected in the intensity spectrum obtained by adding the two signals. A Fourier 
analysis of this spectrum allows for the detection of residual spectral modulation and 
redetermination of ti(A)/t2(A) such that the Fourier power of the modulation dis- 
appears. An alternative method is based on the fact that the RMS of the obtained 
intensity spectrum should be minimal for the correct value of ti(A)/t2(A). 

• The "athermal" multiple-order retarder always has some residual temperature sensitiv- 
ity at least in parts of the wavelength range. These thermal effects, to first order, only 
infiuence the (pL measurement, in a predictable way. However, if the Pl measurements 
are affected by non-ideal effects depending on (pL, the angle needs to be determined 
accurately to apply the calibration to the measurements. A correction can be made 
if the temperature of the modulator during the measurement is known. An indirect 
correction method can be devised in the case that (part of) the measurable 0l(A) is 
constant (e.g. in cases of single scattering). The deviations due to the thermal effects 
can then be easily identified and corrected for over the entire wavelength range. 
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Figure caption listing page 

Fig. 1. (a) Schematic set-up of the spectral modulator. The solid and dashed lines of the 
retarders represent the (orthogonal) fast and slow axis respectively, (b) Illustration of the 
spectral modulation principle on the Poincare sphere. 

Fig. 2. Design of a retrorefiecting quarter-wave retarder prism based on 3 TIRs in a 
fused silica rhomb. The path of light with an angle of 2° from normal incidence is indicated 
in light grey. 

Fig. 3. (a) Thermal test results for the combination of MgF2 (with variable thickness) and 
sapphire as a function of wavelength (550-750 nm). Various thickness ratios (MgF2:Al203) 
from 1.7 to 3.7 have been measured. The vertical axis represents the relative variation of the 
overall retardance of the crystal combination with temperature. The crosses represent the 
determined zero-crossings of the data for a certain wavelength. The dashed line represents 
the thickness ratio 1:2.2. (b) The corresponding theoretical curve with literature values 
for the thermo-optic constants. The dashed line here represents the zero-crossing of the curve. 

Fig. 4. The upper panel shows a modulated spectrum obtained with the prototype 
described in the text (solid curve). The dotted curve represents the unmodulated intensity 
spectrum obtained after adding the signals from two orthogonal orientations of the polarizer. 
The normalized modulation signal is shown below. The lower panel shows the results for 
Pi(A) and (piW as obtained by the retrieval algorithm. Offsets from the ideal values of 
-Pl(A) (=1) and 0l(A) (constant) are clearly observed. The solid curve corresponds to the 
signal shown in the upper panel using the dual beam method, whereas the dotted curve 
shows the results from the single beam method {(pi is the same in both cases). The dashed 
curves represent rotations of the input polarization in steps of 45°, which are reproduced 
by the (pL measurements. The dashed curve for -Pl(A) represents the measurement with the 
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input polarization direction coinciding with the K-prism axis, for which the measurement 
efficiency is significantly higher than for the situation with input at 45° from that (solid 
curves) . 
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